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ABSTRACT 



A force-free surface (FFS) 5 is a sharp boundary separating a void from a region occupied by a charge-separated 
force-free plasma. It is proven here under very general assumptions that there is on 5 a simple relation between the 
charge density fi on the plasma side and the derivative of 5 = E ■ B along B on the vacuum side (with E denoting the 
electric field and B the magnetic field). Combined with the condition 5 — on S, this relation implies that a FFS has 
a general stability property, already conjectured by Michel (1979, ApJ 227, 579): <S turns out to attract charges placed 
on the vacuum side if they are of the same sign as n. In the particular case of a FFS existing in the axisymmetric 
stationary magnetosphere of a "pulsar", the relation is given a most convenient form by using magnetic coordinates, 
and is shown to imply an interesting property of a gap. Also, a simple proof is given of the impossibility of a vacuum 
gap forming in a field B which is either uniform or radial (monopolar). 
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1. Introduction 

An important feature of many models of pulsar mag- 
netosphere is the presence of voids or gaps embedded 
in a charge-separated force-free plasma (HoUoway 1973, 
Ruderman & Sutherland 1975, Rylov 1977, Michel 1979, 
Jackson 1980). Such a void exists in particular in the 
global axisymmetric model first proposed by Michel (1980) 
and later numerically computed by several groups (Krause- 
Polstorff & Michel 1985a, 1985b, Smith et al. 2001, Petri et 
al. 2002, Spitkovsky & Arons 2002). In this latter picture, 
the void extends to infinity, the plasma occupying only a 
bounded region - the electrosphere - around the pulsar. 

In his general discussion of vacuum gaps, Michel (1979) 
reported some explicit examples of such structures, and 
noted that they have all a remarkable property. On the 
surface of discontinuity S separating the plasma of positive 
(negative) charges from the vacuum - a so-called force-free 
surface (FFS) -, the component of the electric field along 
the magnetic field on the vacuum side is always directed in 
such a way that positive (negative) particles be attracted 
by S. Michel made the guess that this phenomenon is gen- 
eral, the existence of a FFS thus implying some form of 
stability. It is the aim of this paper to give a complete ana- 
lytical proof of this conjecture. Actually, we shall establish 
its validity under assumptions much weaker than Michel's, 
as we shall assume neither that the system is stationary 
and axisymmetric, nor that the magnetic field is unaffected 
by the local electric currents. 

The paper is organized as follows. In Sect. 2, we consider 
a general continuous electromagnetic field (E, B) having its 
gradient V(E, B) suffering a discontinuity accross a mov- 
ing surface S, and establish a jump relation for the normal 



derivative of the scalar product E • B. Our starting point 
here is the complete system of Maxwell equations, and a 
set of jump relations for V(E,B) which are established in 
Appendix A. In Sect. 3, the jump relation for 9„(E • B) is 
applied to the case where 5 is a FFS, and it is shown to 
imply a general formula from which the validity of Michel's 
conjecture follows immediately. Sect. 4 is devoted to the 
case of a FFS existing in the axisymmetric stationary mag- 
netosphere of a pulsar. We show in particular how the for- 
mula derived in Sect. 3 can be conveniently rewritten by 
using magnetic coordinates and derive a general property 
of gaps. Incidentally, we give in Appendix B (in addition to 
a simple explicit example illustrating our general relation) 
a proof of the impossibility of a void bounded by a FFS 
and threaded by a field B which is either uniform or radial. 

2. A general discontinuity relation 

Let us consider a domain D of space in which there is a 
moving smooth surface S — S{t). We choose a unit normal 
n on S, and define its "+" side as the one towards which n is 
pointing and its "— " side as the other one. For any function 
X (scalar or vector) continuous in D \S and having well 
defined values on 5^, we define the jump |X| on S to be 



IX] = A^ 



X' 



On S, we introduce the decompositions 

X = X, + X„n 

for any vector field X, and 



(1) 



(2) 



(3) 
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for the nabla operator, with V^ acting parallel to S and 
dn = n • V being the derivative along n. Moreover, in order 
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to define the normal derivatives of X^ and X„ (which are 
yet defined only on S), we extend the vector n - and ac- 
cordingly the decomposition (2) - into a thin layer around 
S by imposing 9„n = n ■ Vfi = (this amounts to set 
fi(r) = n(rp), with Vp the projection of r onto S). Note 
that Vs and 9„ do not commute in general. But obviously 
Vs commutes with |.], i.e., 

lVsXj = V4X]. (4) 

Let us now assume that D contains an electromagnetic 
field (E, B) such that: 

— (E, B) is continuous in D, whence in particular 

[E] = |B1=0 on S. (5) 

— (E, B) is continuously difFerentiable in D \ S, with 
V(E,B) being well defined on S^ and 5+, and 
|V(E, B)] thus being weU defined on S. 

— In D\S, (E, B) does satisfy Maxwell equations 

cV X E = -dtB , (6) 

V • B = , (7) 

V • E = iTT^i, (8) 

cV X B = 47rj -h dtF, , (9) 

where fx and j are the electric charge and current den- 
sities, respectively, and c is the speed of light. 

Under these assumptions, it is possible to derive the 
following set of jump relations for V(E,B): 

[V.El-0, (10) 

[V.B1=0, (11) 

ld„E4 - 47rM , (12) 

|a„B„l - 0, (13) 
An 



[5„B,] 
where 



-7„K[jsl, 



47r 



lllis 



X n, 



7„ 



1 



y2 

* n 



(14) 
(15) 

(16) 



is the squared Lorentz factor associated with the normal 
velocity y„ of S. Moreover, the discontinuity of the normal 
component j„ of the current and that of the charge density 
[I have to satisfy the constraint 

[Jnl = MK. (17) 

Quite certainly, Eqs. (10)-(15) and (17) have already been 
established and used in many other contexts, but we are 
not aware of any references, and then we give a proof of 
them in Appendix A. 

We are now ready to write the general relation express- 
ing the discontinuity of the normal derivative 

|a„(E • B)l = E, • |9„Bsl + S45nS„l 

+B,.|a„E,i + B„|a„£;„i (is) 

of the scalar product E ■ B as a function of Vn and of the 
discontinuities of j and /i. Using equations (12)-(15), we 
obtain 



5„(E.B) 



+—7' n X E 



Vn 



B, 



•iJs 



(19) 



3. Application to a force-free surface 

3.1. Assumptions 

We now assume that 5 is a FFS, i.e., there is a charge- 
separated force- free plasma on one side of 5 - S~ , say - and 
a vacuum on the other side. We thus have the additionnal 
bulk relations 



V X B 
E = 



inside the plasma, and 
j = 0, 

^I = 0, 

inside the vacuum. This implies 

[j] = -M^v- , 

Eg = -Bn — X n n x B^ 

c c 



(20) 
(21) 



(22) 
(23) 



(24) 
(25) 

(26) 



on iS, where we have noted to get the last relation that, as 
a consequence of Eqs. (17) and (24)-(25), i;,7 = U„ - this 
equality just expressing the fact that S stays a plasma- 
vacuum interface at any time. 

3.2. The discontinuity relation 

Hereafter, we denote as i an arclength along a magnetic 
line, with i increasing when passing from the plasma to the 
vacuum, and as u the associated tangent unit vector (then 
u • n > on tS). Our aim is to compute the value of the 
derivative with respect to (. of the component 



E 



E u 



of the electric field. 

As E • B = on the plasma side, we have on S 

E-B = 0, 

[V,(E.B)]+ = 0, 
ia„(E.B)l = [a„(E.B)]+. 

Then E\\ satisfies on S 



E^ 



S|i=0, 



(27) 



(28) 
(29) 
(30) 



(31) 



and 



B 



2<^ 

d£ 



= [B-V(E-B)] + 

= [B, •V,(E.B)] + 
= B„I5„(E.B)1, 



[B„a„(E . B)] 



(32) 



where we have used equations (29)- (30) to get the last 
equality. Using Eqs. (24)- (26) in the right-hand side of equa- 
tion (19) and remembering that v^ — Vn, we get 



ia„(E.B) 



-Att^ B„ (1 - 7, 

-A'Kfi^Bnlll^'^ , 



:{Vs)' 



(33) 
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where 



7 =1 ji 



(34) 



By using the latter relation in equation (32) , we obtain 
finally 






R2 

. _ ^n 2 -2 

= "47rM -g27„7 : 



or, by introducing the angle a between B and n, 



dEi\ 

~de' 



= —Anji cos^ 0^7^ 7 ^ ■ 



(35) 



(36) 



On a part of S not tangent to B (cos a ^ 0), we thus 
have 



~dl 



= — 47r(/i cosQ!7„7 ) < 



(37) 



which implies along with Eq. (31) the existence of a layer 
on the vacuum side in which 



H'Eu <0, 



(38) 



i.e., the sign of E'n is opposite to that of fj, . When /x > 
0, E\\ is directed towards S, which attracts the positive 
particles, while when /i^ < 0, E^\ is directed away from 5, 
which attracts the negative ones. This is just the general 
stability property conjectured by Michel (1979). 

As an illustration of our general formula (36), we con- 
sider in Appendix B a particular example of FFS borrowed 
from Michel 1989. 



4. Application to the aligned pulsar magnetosphere 

4.1. Assumptions and the basic discontinuity relation 

In this section, we use standard spherical coordinates 
{r,9,(p) attached to the Cartesian frame (0,x, y,z). We 
consider a pulsar of center O rotating at the angular veloc- 
ity i7,z and surrounded by a stationnary magnetosphere D. 
The latter is axisymmetric around the z-axis - this allows 
us to consider the equations governing its structure as be- 
ing set in the meridional half-plane 11 = {(/? = 0} -, and it is 
made of two regions separated by a FFS S: Dp, which con- 
tains a charge-separated force-free plasma, and Dy which 
is a vacuum. Here we allow D^ and D^ to be constituted 
of several connected parts and/or to extend to infinity (we 
do not refer to any specific model). We denote as F and dT 
the intersections of D^ and S, respectively, with H. 

In the whole D: (i) the magnetic field admits the well 
known representation 



B 



V^ X if 



(39) 



with the flux function A satisfying 

An 



in Dy (with /x = 0) and in the part of Dg where the parti- 
cles do not stream along the lines and then have a purely 



rotational motion at the angular velocity 17; (ii) the electric 
field can be written in the form 



E = -VV' 



(41) 



owing to the stationarity assumption, with the electrostatic 
potential ip solving Poisson equation 

- VV = 47r^ • (42) 

Inside the plasma, ip — ipiA), with 



^AiA) = 



dV- _ O 



(43) 



(see, e.g., Goldreich & Juhan 1969). 

Let us write the particular form taken by Eq. (36) on 
dr. As S does not move: (i) l^i — and 7„ = 1; (ii) a 
charged element on S has a pure rotational motion at the 
angular velocity fl given by Eq. (43) {ft — il, if the element 
is connected to the star by a piece of magnetic line fully 
embedded in Dp). We thus have 



d^ll 
'dl 



d^V' 
d£2" 



-4:11 fi 7 ^ cos^ a , 



with 



7^2 = 1- {rsmOn/cY' 



(44) 



(45) 



4.2. Expression of the discontinuity relations in magnetic 
coordinates 

Inside the vacuum region Dy, the magnetic field also admits 
the representation 



B = Vy, 

with the scalar potential V satisfying 
VV = 0. 



(46) 



(47) 



It is often convenient to use A and V as curvilinear orthog- 
onal coordinates in F (note that VA • VV^ = 0). The vac- 
uum electrostatic potential written as a function "0(^4, V) 
thus satisfies 



V?A = ij'A^A + ij'y'^V , 



E-B = -^'vB 



(40) V^^ 



2 

^ [(rsin0)VAA + ^yy]S'+2V'AS. = 

(Eq. (40) has been used to get the last line). 
By Eqs. (49) and (28), we have on dT 

V-y-O, 

while Eqs. (49) and (44) imply 

47r^"7"2 cos^ aB'^ = -B • [V(E • B)] + 
= B . [V{B^i;'y)\+ -- 

whence 

47r/i^7^^ cos^ a 



i?V^y 



vv 



B2 



(48) 
(49) 

(50) 



(51) 



(52) 



(53) 



It is interesting to note that the two other second deriva- 
tives of V'(^: ^) can be also given explicit expressions. To 
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get V'ai/' ^^ differentiate Eq. (51) along dT. Denoting by I 
an arclength increasing in the direction of the tangent unit 
vector t ~ X h and vi^ith a being oriented by ip, we have 



dA 



= ^'Avr-n + ^vv\ 



dl 



dv_ 



^AV I ^ ^^"^ ^ '^'^^ '^^ + ^VV I ^^"^ '^^ ■ 



and then 



^v 



47r/i 7 ^ sin a cos a 
r sin^B^ 



(54) 



(55) 



And by using Eq. (50), we obtain 



^'; 1+ 



-1 



''AA 



(rsin6lS)2 

-47r 



(rsin 



(47r^"7~2 cos^ a + 2%Ij'j^B^) 
^"7"^cos2a-^,77^— ) , (56) 



where /i* = —VI^:Bz^1/2ttc is the Goldreich- Julian or coro- 
tational charge density (Goldreich & Julian 1969) and 7* is 
given by Eq. (45) with fl = fi*. 

Near the pulsar, there is a region D' where the particles 
move at a velocity small compared to c and the magnetic 
field is essentially determined by the sole currents flowing 
inside the star (therefore both A and V are a priori known 
quantities in I?', and using them as coordinates therein 
seems to be especially appropriate). In many papers, the 
authors consider a FFS which is located in such a region 
(this is the case in all the calculations of Michel's electro- 
sphere referred to in the introduction) and it is possible to 
take 7 = 1 in all the relations above. In that situation, a 
simplified proof of Eqs. (53), (55) and (56) can be obtained 
by making the nonrclativistic approximation from the very 
beginning. In fact, it is just needed to use Eqs. (40) and 
(42) (the former with /i = everywhere) and to take into 
account the continuity of tjj and VV' on S. 

4.3. A general property of gaps 

Let us define a vacuum force-free surface (VFFS) to be a 
surface located inside a void on which E • B = 0. Then 
we have the following result: there is necessarily such a 
surface in _D„ if the latter contains a piece of a magnetic 
line connecting two points ri and r2 of dT at which the 
charge density has the same sign - i.e., /i(ri)/i(r2) > - 
and cos a 7^ 0. 

The proof of this statement is immediate. As a conse- 
quence of Eqs. (28) and (36), we have 



,5(ri) = ,5(r2) = 0, 

where (5 = E • B, and 

(B-V(5)(ri).(B.V<5)(r2)<0, 



(57) 



(58) 



which clearly implies that 6 has to vanish at least once 
(in fact 2p + 1 times, with p > an integer) on the line. 
By continuity, the same result holds on all the neighboring 
lines, and the points where S = form a surface indeed - a 
VFFS. 

A trivial illustration of this result is provided by the 
calculations of Michel's electrospheric structure in a purely 
dipolar field. In that case, there is a bundle of magnetic 
lines connecting through the unbounded vacuum region the 



negatively charged domes overlying the north and south 
poles, respectively, and it is a priori quite obvious that this 
bundle threads a VFFS which is just the part of the equa- 
torial plane it does intersect - we have indeed E • B = 
on{z = 0}nDasa consequence of the symmetry of the 
system. 

As a more interesting illustration, we may consider - 
following Asseo et al. (1984) - the theoretical possibility of 
a connected gap dy fully included in the negatively charged 
region overlying the North pole, say. If such a gap does 
exist, it necessarily contains by the result above a VFFS 
Sy (possibly made of several parts) cutting all the magnetic 
lines an odd number of times. In particular, Sy meets the 
2;-axis at a critical point of "0 (VfA — 0) where at least 
two equipotential surfaces meet (this point has to be of 
the X-type as the vacuum potential cannot reach either 
a maximum or a minimum inside dy). As E also vanishes 
at the points of intersection of the boundary of dy with 
the z-axis, with each of these points generating a cusped 
equipotential surface (Asseo et al. 1984), we see that the 
potential needs to have a quite curious structure, which 
makes the existence of dy unlikely (this point is discussed 
in more details in a forthcoming paper, where we combine 
the result here with other properties of the potential ip to 
recover Asseo et al.'s nonexistence result without appealing 
to physical arguments involving pair production). 

Incidentally we note that the impossibility of a station- 
ary vacuum gap (either of finite extent, or unbounded with 
E vanishing at infinity) can be most easily proved in the 
nonrclativistic approximation when the given potential field 
B is uniform or radial (see Appendix B). An hypothetical 
pulsar with a monopolar (or "split monopolar") magnetic 
field could then not have a finite electrosphere of Michel's 
type. It is worth noticing that Michel (1973) has found in 
that situation an exact gapless solution of the Goldreich- 
Julian type, with a wind extending to infinity. 

5. Conclusion 

In this paper, we have provided a general analytical proof 
of a stability conjecture put forward by Michel (1979) some 
twenty five years ago. According to the latter, any FFS does 
attract particles located on the vacuum side if they are of 
the same sign as the particles on the plasma side, a property 
which makes the electrostatic at work in a pulsar magneto- 
sphere with gaps somewhat nonstandard. Although Michel 
formulated his conjecture on the basis of particular axisym- 
nietric, stationary and nonrclativistic examples, we have 
shown here that it was quite robust, still holding true when 
neither of these simplifying assumptions are made. 

In the case where the FFS exists in the magnetosphere 
of an "aligned pulsar" , we have shown that the general re- 
lation from which we have derived the validity of Michel's 
conjecture takes a particularly simple form when magnetic 
coordinates are used. This relation determines indeed the 
second derivative of the vacuum electrostatic potential ■0 
with respect to the magnetic potential \^ as a function of 
the charge density, the condition that 5 be a FFS being 
expressed by the vanishing of the first derivative of ip with 
respect to that same variable. We have also derived an in- 
teresting property of a vacuum gap threaded by magnetic 
lines connecting plasma elements having the same sign of fi 
- such a gap does necessarily contain a VFFS - and proved 
the impossibility to have in a given uniform or radial mag- 
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netic field a bounded stationary gap, or an unbounded one 
with no electric sources at infinity. 

Appendix A: Jump relations for V(E,B) 

We give here a derivation of the jump relations for V(E, B) 
quoted in Section 2. 

1. Derivation of Eqs. (10)- (11). As an immediate conse- 
quence of Eqs. (4)-(5), we have 



where we have taken Eq. (A. 7) into account, and 

9„EJ . (A.17) 



n X 



^B.] = ^[W-^1 



Finally, Eqs. (A. 15) and (A.17) combine to give 

Air 



5„E,1 



|V,E], IV.B]) = (V4E], V4B]) = (0,0) 



(A.l) 



6 



^llVnM 



(A.18) 
(A.19) 



2. Derivation of Eqs. (12)-(13). By the standard decom- 
position 



with 7„ defined by Eq. (16). 



V • X = V, • X + a„x„ 



(A. 2) Appendix B: Particular configurations 



of the divergence operator on a surface (e.g.. Brand In this appendix, we present a brief study of three particular 

configurations. We assume stationarity - then E can be 
expressed in terms of an electrostatic potential -0 - and 
nonrelativistic velocities, and take B to be a given potential 
(A. 3) field. We use Cartesian coordinates (x, y, z) and associated 
spherical coordinates {r,0,(p). 



1947) and Eq. (A.l), we have 
IV-(E,B)1 = lVs-{E,B)] + ldn{En,B„ 



On the other hand, we have by Eqs. (7) and (8) 

IV-Bl = 0, 
[V • El = AttM , 

and then 

[a„B„i = 0, 

ldnE4 - 47r|Ml . 



(A.4) 
(A.5) 



B.l. Spherical plasma confined by a uniform B 



3. Derivation of Eqs. (14)-(15) and (17). We note that Eq 
(5) implies 



Following Michel 1989, we consider a ball V of center O 
and radius tq containing free electric charges distributed 
with the uniform density fi. V is threaded by a uniform 
magnetic field Bz [B > 0) and it is submitted to an external 
(A. 7) quadrupolar electric field oriented in such a way that the 
total electric field be given by 



(A.6) 



iTTjir^ 



d4E,Bl = ia,(E,B)l +K[a„(E,B)] =0, 



(A.8) E : 



where dj is the time derivative following the normal mo- 
tion of iS, and we use the decomposition 



„-5-^r + /c(xx + yy — 2zz) for r > tq , 



^^r -I- fc(xx -\- yy — 2zz) for r < tq . 



(B.l) 



VxX = VsXX + nx dn^s 



(A.9) 



of the rotational operator on a surface (Brand 1947) and 
Eq. (A.l) to write 

IV X (E, B)l = |V, X (E, B)] + n X ia„(E„ B,)] 

- (nxia„E,l,nx|a„B,]). (A.IO) 

Plugging these relations into the ones 

c|VxEl = -I9,B1, (A.ll) 

c|V X Bl = 47r|j] + IdtEj , (A.12) 

resulting from Eqs. (6) and (9) thus leads to 

cnx|a„E,i = v;ia„B], (A.13) 

en X ia„B,l = 47r|j] - K[a„E] . (A.14) 

The normal component of equation (A.13) gives once 
more Eq. (A.6), while its tangential component is equiv- 
alent to 



From now on, we choose the parameter k = 2Tr fi/3. Then 
E • B = BEz = in F, and the plasma turns out to be 
force-free if we take it to move at the drift velocity 



Ex B 



27Tcr sin On _ 
B ^ 



(B.2) 



corresponding to a rigid rotation at the angular velocity 
fi = —2ttcij,/B. Of course, E^ — on the sphere {r = rg} 
separating the plasma from the vacuum - it is a FFS indeed. 
In the vacuum, 



[a„E, 



Vn 



IdnBs 



(A.15) 



and 

dE, ^ 
dz 

whence 
dE 



Attjiz f rQ 






Anfi 
~3~ 



- 1- 



irf^z^ 



By taking similarly the normal and tangential compo- 
nent of Eqs. (A.14), we obtain 



dz 



(ro , 0) = —Anfi cos = — 47r/i cos a : 



(B.3) 



(B.4) 



(B.5) 



[jnl = MVn 



(A.16) 



we recover here in a particular situation our general relation 
(36) with 7 = 7„ = 1. 
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B.2. Nonexistence of voids in some particular B 

Let D^ be a void in contact with a force-free plasma region 
De all along its boundary S - a. FFS. D^ may be either of 
finite extent, or unbounded, in which case we also require 
that the electric field vanishes at infinity. 

Wc first consider the case where B is uniform (B = Bz). 
Then 



V^E, = 


in Dy , 


(B.6) 


E,^0 


on S , 


(B.7) 


lim ^2=0 


(if Dy unbounded) , 


(B.8) 



from which we can conclude at once by a standard unique- 
ness theorem for Laplace equation that E^ = —dzip = 
in Dy. Therefore, ip = f{x,y), with / an harmonic func- 
tion defined on the projection of Dy onto the plane (x, y). 
Because the potential is continuous on S and is constant 
along the lines in Dg, we must have ip{x,y,z) — f{x,y) 
and then /x = — V^//47r = in that part of D^ magneti- 
cally connected to Dy. This is in clear contradiction with 
the definition of D^, and we can conclude that the gap Dy 
cannot exist. (Note that there is no disagreement between 
this result and the example of the previous subsection, as 
in the latter Ez is not bounded when r — ;> oo). 

Next we consider the case where B is radial (B is the 
monopolar field B = Bo{ro/r)'^f, or a split monopolar field 
with zero flux, but with a current sheet) and Dy is located in 
the region r > ro, say. We note that the function x = r-'Vip 
satisfies 



v\ = o 


in Dy , 


(B.9) 


X-0 


on S , 


(B.IO) 


lim X - 


(if Dy unbounded) , 


(B.ll) 



which implies that x = rdr'>p = in Dy. Then tp — g{9, ip) in 
Dy, with g harmonic, and we can conclude by an argument 
similar to the one above that this relation also holds true 
in a part of D,, - a contradiction. Vacuum gaps do not exist 
in the radial field either. 
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